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Abstract. We analyze the resonant inelastic x-ray scattering (RIXS) spectra at 
the Cu and Ni K edges in La2Cu04 and La2Ni04, respectively. We make use of 
the Keldysh- Green-function formalism, in which the RIXS intensity is described 
by a product of incident-photon-dependent factor and density-density correlation 
function in the 3d states. The former factor is calculated using the 4p density of 
states given by an ab initio band structure calculation and the latter using the 
wavefunctions given by a Hartree-Fock calculation of a multi-orbital tight-binding 
model. The initial state is described within the Hartree-Fock approximation and 
the electron correlations on charge excitations are treated within the random phase 
approximation. The calculated RIXS spectra well reproduce several characteristic 
features in the experiments. Although several groups have interpreted the RIXS 
peaks as bound excitons, our calculation indicates that they should be interpreted 
as band-to-band excitations augmented by electron correlations. The difference 
in RIXS spectra between La2Cu04 and La2Ni04 is explained from this point of 
view. 
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1. Introduction 

Investigation of the electronic excitations in the transition metal compounds is 
fundamental for understanding their electronic properties. The excitations may be 
characterized into two types, spin and charge excitations. For the former, the inelastic 
neutron scattering technique is quite useful to investigate energy- momentum relations. 
For the latter, taking advantage of strong synchrotron sources, the resonant inelastic x- 
ray scattering (RIXS) technique has become a powerful tool to probe charge excitations 
and to investigate their energy-momentum relations. For the RIXS tuned at transition 
metal if -edge, the process is described as a second-order optical process, in which a 
Is-core electron is excited to an empty 4p state by absorbing a photon, then charge 
excitations are created in the 3d states to screen the core-hole potential, and finally 
the photoexcited Ap electron recombines with the Is-core hole by emitting a photon. 
In the end, the charge excitations are left behind. 

RIXS measurement at the Cu K-eAge has been carried out to search for charge 
excitation modes and to determine the momentum dependence of charge excitation 
in high Tq cuprates and their related materials for better understanding of the 
unconventional high-Tc superconductivity. For the undoped material La2Cu04, the 
RIXS spectra as a function of energy loss are composed of several peaks and show 
clear momentum dependence[ll [21 [31 [4]. The peak position in terms of energy loss 
shifts remarkably as the momentum transfer changes. Several groups have claimed 
that the RIXS intensity refiects a bound exciton and interpreted the peak-shifts as 
a dispersion relation of bound excitonfH [^. The RIXS experiment at Ni K-eAge on 
La2Ni04, which has isostructure with La2Cu04, has been also carried out in order to 
clarify the difference of charge excitations from La2Cu04. In contrast to the cuprate, 
no superconductivity has been reported in the nickelate. It has been reported that 
the RIXS spectra are composed of several peaks which do not show clear momentum 
dependence. The peak position in terms of energy loss does not shift even though the 
momentum transfer varies. It has been suggested that the exciton is well localized in 
the nickelate, while it is mobile in the cuprate[3j. 

On the other hand, it has been argued that the RIXS intensity should be 
interpreted as band-to-band excitations augmented by electron correlations, not as 
bound excitons. Nomura and Igarashi (NI)(6]|7j have proposed a general formalism 
of the RIXS spectra by extending the resonant Raman theory developed by Nozieres 
and Abrahams[5 on the basis of the many-body formaHsm of Keldysh. In the NI 
formula, with the assistance of the Born approximation on the interaction between 
the Is hole and 3d electrons, the RIXS intensity is described by a product of an 
incident-photon-dependent factor and a density-density correlation function in the 3c? 
states. The NI formula has been successfully applied to the quasi-one-dimensional 
cuprates SrCuOs^fij, CuGeOalQj, two-dimensional cuprate La2Cu04[71 [ini [IJ, and 
the prototypical AFM insulator NiO[l2]. In these studies, the electronic structures in 
the AFM phase have been calculated within the Hartree-Fock approximation (HFA) 
using tight-binding models. It is known that the HFA works well for describing 
the electronic structures in the AFM insulators. Two-particle correlations have been 
taken into account within the random phase approximation (RPA), which has been 
found to modify the spectral shape as a function of energy loss, having led to a good 
agreement with the experiments. On the basis of these successes, we may conclude 
that the RIXS intensity arises from band-to-band transitions to screen the core-hole 
potential in the intermediate state. Multiple-scattering contributions due to the core- 
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hole potential have been also investigated in order to examine the validity of the 
Born approximation, because the core-hole potential is not definitely weak. Having 
evaluated the contributions by means of the time-representation method by Nozieres 
and De Dominicis[i3j, it was found that the contributions could be mainly absorbed 
into the shift of the core-level energy with minor modifications of the RIXS spectral 
shape[T^. This result partly justifies the use of the Born approximation. 

In the present study, we discuss the RIXS intensity and clarify the origin of the 
difference of them in La2Cu04 and La2Ni04 by applying the NI formula. We exploit 
the tight-binding HFA calculation and the ab initio band structure calculation based 
on the local density approximation (LD A) to calculate the density-density correlation 
function and the incident photon dependent factor, respectively. We obtain stable 
AFM insulating solution having an energy gap aboutl.7(3.5) eV within the HFA 
for La2Cu04 (La2Ni04). Note that the band structure calculation with the LDA 
fails to describe the AFM insulating state for these systems. The two-particle 
correlations in the intermediate states are treated within the RPA. We obtain the 
RIXS intensities in a range of energy loss 2 6 (4 '--^ 8) eV for La2Cu04 (La2Ni04). 
The calculated RIXS spectra well reproduce several characteristic features in the 
experiments [2131 H]. Particularly, the calculated RIXS spectra for La2Cu04 exhibit 
the noticeable momentum dependence of the peak structure while those for La2Ni04 
do only slightly, being consistent with the experiments. Although our results do not 
reproduce all features in the RIXS spectra at this stage partly because of the strong 
correlation between 3d electrons, our results suggest that the RIXS intensity should be 
interpreted as band-to-band transition and the peak shifts as a function of momentum 
should not be interpreted as a dispersion relation of a kind of exciton, because the 
peaks with broad width constitute an energy continuum generated by band-to-band 
transition. 

The present paper is organized as follows. In Sec. 2, we briefiy summarize the 
NI formula for the RIXS spectra. In Sec. 3 we discuss the calculated RIXS spectra 
comparing with the experiments in connection with the electronic structure within 
the HFA in the AFM phase of La2Cu04 and La2Ni04. The last section is devoted to 
the concluding remarks. 

2. RIXS Formula 

We briefiy summarize the NI formula for the RIXS. In the RIXS process, an incident 
photon is absorbed by exciting a Cu(Ni) Is core electron to the unoccupied Cu(Ni) Ap 
state, and a photon is emitted by recombining the 4p electron and the core hole. In the 
intermediate state of the RIXS process, the core-hole potential acts on the 3d states 
and creates electron-hole pairs in order to screen the core-hole potential. In the end 
of the process, single electron-hole pair is left behind within the Born approximation, 
carrying momentum-energy q = (q, Ci^) = (q^ — q/,^;^ ~ujf), where qi — (q^jW^) and 
9/ = momentum-energy of incident and scattered photons, respectively. 

The RIXS intensity is derived on the basis of the Keldysh-Green function scheme. 
The process is diagrammatically shown in FigUl Within the Born approximation to 
the core-hole potential, the RIXS intensity is given by. 




{q)JBXX' {ujt,e,;ujf,ef) , 
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Figure 1. Diagrammatic representation for ttie RIXS intensity. The wavy and 
solid lines represent photon and electron bare Green's functions, respectively. The 
dotted line is the core-hole interaction V . The solid lines with double arrows are 
the Keldysh type Green's functions. The shaded area represents the effective 
scattering vertex renormalized by 3a!-3(i Coulomb interaction in the RPA. 



where A indicates the Cu(Ni) site A in a unit cell. TV is the number of unit cell and w 
is the ls-4p dipole transition matrix element, which is assumed to be constant. The 
factor Jbaa' (^i, e^; w/, e/) describes the incident-photon dependence, which is given 

by 



fn' 



riri' 



X 



E^"'^BA' , (2) 



where e^^ (e/?;) is the ry component of the polarization vector (e/) with rj = x,y, z, 
-m' 



and Lg*^ (LUi;ujf) is given by 



V represents the Is core-hole potential acting on the id electrons, Tis is the life-time 
broadening width of the core-hole state, and the lower limit of the integral eo indicates 
the energy at the bottom of the conduction band. This expression comes from the 
upper triangle in Figdl The p^^^ is the DOS matrix in the p symmetric states at the 
Cu(Ni) site A, which may be given by 

P% i^)-T.T. '^A.. in, k) cl^wa {n, k)6{e- e„ (k)) , (4) 

(T nk 

where ^aijo- {n, k) is the amplitude of component with spin a at the Cu(Ni) site A 
in the band state specified by the band index n and momentum k with energy e„ (k). 
The factor i^^rn'o-' Xma il) ^1- O density-density correlation function of the 

Keldysh type, which is defined by 

/oo 
((PqA'™'.')^(r)PqA™.(0))e-MT, (5) 
-oo 

where 



PqAmj — \J j^^j ^k+aAmcr^kAin<T, (6) 
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dk\mc7 is an annihilation operator for the state denoted by the index kAmcr. The index 
Xma specifies a tight-binding orbital at the Cu(Ni) site A with 3d orbital m and spin 
a. Wavevector k in eq. ^ runs over the first Brillouin zone. Vectors r„ and ua in 
eq. |[7]) represent a position vector of the n-th unit cell and the position vectors of the 
Cu(Ni) site A in the unit cell n. 

We calculate the density-density correlation function ^ with taking account of 
the correlation effect on the electron-hole pair by the RPA as shown in fig. [TJ For 
more details of the derivation, see refs. [TO t fT2l ITT] . 

3. Results and Discussion 

In order to calculate the incident-photon-dependent factor Jbaa' (t^ij e^; w/, e/), we 
need the 4p DOS pj^p{e) on the Cu(Ni) site A. We use the 4p DOS given by the band 
structure calculation based on the LDA assuming non magnetic states. Although the 
LDA fails to describe the antiferromagnetic insulating states with a wide energy gap, 
the calculated Ap DOS in the conduction band may be reliable because the 4p DOS 
has large intensity in the energy range above about 10 eV higher than the top of 
the occupied states. Due to the crystal symmetry, pj^pie) is almost independent of 

Cu(Ni) site A, P\4p{e) ~ P\ip(^) P\4pi^)' ^^'^ the values for r\ ^ r( are negligible. 
We set the energy difference between the Cu(Ni) Is level and the prominent peak 
in Pxipi^) to be 8992 (8348) eV as shown in fig. [H Under the condition that the 
dipole matrix element is constant and that the interaction between the core hole and 
the 4p electron is neglected, the 4p DOS becomes proportional to the Cu(Ni) K-edge 
absorption spectra. In the following, we assume that the incident and emitted x-rays 
are polarized along the z direction. Accordingly, the factor Jbaa' ('j-'i: e^; w/, ey) is 
replaced with the factor \L^\ {uji;uj f )\'^ . 

Another factor Yy^'a Xma(l) mainly determines the structure of the RIXS 
spectra as a function of energy loss. We calculate this factor using the tight-binding 
wavefunction given by the HFA assuming the AFM order. The Cu(Ni) 3d orbitals, 
apical and in-plane O 2p orbitals are included in the tight-binding model. The 
Slater-Koster parameters are taken from the LAPW band structure calculation for 
La2Cu04[l4]. The same values are assumed for La2Ni04, because reliable Slater- 
Koster parameters are not available at this stage. Although a precise calculation 
based on the first principle is strongly required, this assumption is plausible because 
the bond lengths Cu — O and Ni — O are nearly equal with each other[l5l [16] and the 
band structures near the Fermi level calculated within the LDA assuming non magnetic 
states for La2Cu04 and La2Ni04 are similar to each other|17|. The intra-atomic 
Coulomb interaction on the Cu(Ni) sites is given by the Slater integrals F^, and 
F'^. Among them, and F'^, which are known to be slightly screened by solid-state 
effects, are taken from the analysis of the x-ray photoemission spectroscopy ^Hl HI]. On 
the other hand, F^ is known to be considerably screened, so that we regard the value 
as an adjustable parameter. The Coulomb interaction on O sites and that between 
Cu(Ni) 3d and O 2p orbitals are absorbed into a renormahzation of the O 2p level 
parameters. The Cu(Ni) d-level position relative to the O p-levels is given by the 
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Figure 2. pj^p and + /'lilp j comparing with the 

fluorescence yield (gray solid curve) reproduced from ref. (1]. Energy zero of the 
DOS is shifted so that the energy of the prominent peak in p^^ coincide with the 
fluorescence peak. 



charge-transfer energy A defined as A = — -Ep + 9L/ {IS. — Ed — Ep + 8U) in the 
dp {(fi) configuration for La2Cu04 (La2Ni04), where U is the multiplet-averaged d-d 
Coulomb interaction given hy U = F° - (2/63) F^ - (2/63) F'^[2Qj. The charge-transfer 
energy A is also treated as an adjustable parameter in the present calculation. We 
assume C/ = 11 eV, A = 2.5 eV (C/ = 6 eV, A = 4 eV) for La2Cu04 (La2Ni04). 

Figure [3] shows the DOS projected on the Cu(Ni) — r^ and — y^ states. The 
calculated energy gap is about 1.7 (3.5) eV for La2Cu04 (La2Ni04). The — y^ states 
of the local minority spin reside on both occupied and unoccupied energy regions with 
noticeable weight. Therefore, the transitions from the — states in the energy 
ranges denoted by B, C, and D to the x^ — y^ states in the energy range denoted by 
A contribute to the screening process of the Is core-hole potential, consequently to 
the RIXS intensity. We note that the transitions between the states having different 
local symmetry from each other is forbidden. The distribution of the x^ — states of 
the local minority spin is similar between La2Cu04 and La2Ni04 except for the band 
width. The band width of the states including x^ — y^ state of the local minority spin 
for La2Ni04 is much narrower than that of La2Cu04 due to the larger charge-transfer 
energy A. In La2Cu04, the Zz^ — r"^ states are completely occupied for both spins . In 
La2Ni04, a small amount of the 3z^ — states of the local minority spin reside even 
in the occupied energy region, while those of the local majority spin are completely 
occupied. Thus, the transition in 3z^ — channel dose not occur in La2Cu04, while 
the transition can contribute to the RIXS intensity in La2Ni04. In contrast to the eg 
states, the tig states cannot contribute to the RIXS intensity, since the local majority 
and minority spin tig states are completely occupied. 

Combining the two factors, we obtain the RIXS spectra with assuming the core- 
hole lifetime broadening Fi^ — 1 eV. We convolute the result with a Lorentzian 
function of FWHM = 0.3 eV for taking account of the instrumental resolution. Figure 
[Ushows the calculated spectra for La2Cu04 at several momentum transfer q comparing 



Charge Excitations in Cuprate and Nickelate in Resonant Inelastic X-Ray Scattering? 



0.8 



0.8 



- Cu 


3z -r 

x'-y" _ 

A 


- IL 

d| c b 


1 





0.8 


- Ni 




1 ' 











A 








f 


















0.8 


- D C B 


1 







I I \ I I I L 

-4-2 2 

Energy (eV) 



Figure 3. Partial DOS of — and — orbital at a transition metal 
site and of pa orbital at an in- plane oxygen site for La2Cu04 (upper panel) and 
La2Ni04 (lower panel). Energy zero is at the top of valence band. DOS for the 
local majority and the local minority spins are presented in upper and lower half, 
respectively. Roman letters A, B, C, and D are denoted to the energy ranges 
where the — partial DOS in the local minority spin has noticeable intensity. 



with experiments[3]. The RPA correction modifies the spectra calculated within HFA, 
leading to better agreement with the experiments. We obtain continuous spectra as a 
function of energy loss to ranging from w = 2 eV to 6 eV. Intensities around lu — 2.2, 
3.2, and 4.5 eV are caused by charge excitations of B^A, C— >A, and D^A transitions, 
respectively, within the x"^ — y"^ symmetry in the local minority spin states (see fig. 
[HI). A prominent peak around 4.5 eV stays at the same position with changing the 
momentum transfer q. The spectral shape in the low energy region changes as the 
momentum transfer q varies. A broad hump existing around lu — 2 — 4 eV at q = (0, 0), 
which is enhanced by the RPA correction, grows up to become a peak around 3.2 eV 
at q = (7r,0). This behavior corresponds well to the experimental fine shape at 
q = (tt, 0), which looks like two peaks around uj — 3.2 and 4.5 eV. On the other 
hand, the intensity of the broad hump around (jj = 2 4 eV is suppressed by the 
RPA correction at q = (tt, tt) and only one peak is overwhelmingly left at 4.5 eV. This 
explains the experimental spectra at q = (tt, tt), which looks like a single peak around 
u) = 4.5 eV. Although several characteristic features of RIXS spectra are reproduced in 
a good agreement with the experiments|4j|3j, several discrepancies between the present 
calculation and the experiments still remain; the 4.5 eV peak is hardly discernible at 
q = (0,0), the intensity at w = 4 eV is enough large forming a peak-like structure at 
q = (0, 0) with = 8992 eV[2]. To remove these discrepancies, we may need to take 
account of the electron correlations beyond the RPA and the effects beyond the Born 
approximation to the core-hole potential. 

Figure [H shows the calculated spectra for La2Ni04 at several momentum transfer 
comparing with experiment [J- We obtain continuous spectra ranging from = 4 
eV to 8 eV. The RPA correction works to suppress the spectra calculated within the 
HFA in a energy range 4 ~ 7 eV at all momentum transfer q. The intensities around 
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Figure 4. Calculated RIXS spectra (left) and experimental spectra (right) from 
ref. for La2Cu04. The experimental spectra are reproduced from ref. Solid 
and dotted curves are .the spectra calculated with RPA and HFA, respectively. 
Incident photon energy is assumed to be 8992 eV, although it is 8992.5 eV in the 
experiment. 




Figure 5. RIXS spectra for La2Ni04. The experimental spectra are reproduced 
from ref. [1]. Incident photon energy is assumed to be 8348 eV. 



w = 4, 5, and 6.3 eV arise mainly from the charge excitations of B^A, C^A, and 
D^A transitions, respectively, within the — symmetry in the local minority 
spin states. The transition between the — states of the local minority spin also 
contributes to the intensity around u = 5.8 and 7.3 eV. In contrast to La2Cu04, the 
calculated RIXS spectra show only sHghtly peak-shifts and intensity-changes as the 
momentum transfer q changes. This non-dispersive behavior of the peak-structure as 
a function of energy loss is well consistent with the experiment [Ij. We note that this 
is a consequence of the fact that the band width of the — states of the local 
minority spin in La2Ni04 is much narrower than that in La2Cu04 due to the larger 
charge-transfer energy A. 
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4. Concluding Remarks 

We have analyzed the momentum transfer dependence of the RIXS spectra in La2Cu04 
and La2Ni04 on the basis of the formula developed by Nomura and Igarashi. This 
formula expresses the RIXS spectra by a product of the density-density correlation 
function and the incident-photon dependent factor. This formula makes it possible to 
calculate the RIXS spectra on the complicated models including many orbitals, and 
provide clear physical interpretations to the RIXS spectra. It also may be possible to 
incorporate this formula into an ab initio band structure calculation up to the RPA 
level, because the formula is similar to that of the dielectric function e (q, uj). In the 
present study, a multiorbital tight-binding model, which includes all the Cu(Ni) 3d 
and 2p orbitals as well as the full Coulomb interaction between 3d orbitals have been 
used to calculate the density-density correlation function, which has been calculated 
from the tight-binding wavefunctions within the HFA and RPA. The incident-photon- 
dependent factor has been evaluated from the 4p DOS given by the ab initio band 
structure calculation. 

Our results reproduced several characteristic features in the RIXS experiments as 
a consequence of band-to-band transition between the occupied and unoccupied x^ — y^ 
states of the local minority spin. Present calculation has also shown that in La2Ni04 
the RIXS peak position in terms of energy loss shifts only slightly as the momentum 
transfer varies, while that looks to shift remarkably in La2Cu04, being consistent 
with the experiments. This is simply explained as a consequence of the fact that the 
band width of the states contributing to the RIXS intensity is significantly narrower 
in La2Ni04 than in La2Cu04 because of a large charge-transfer energy A. Present 
calculation does not reproduce all features at this stage partly because of the strong 
correlation between 3d electrons. In spite of several drawbacks of our analysis, it is 
suggested that the peak shifts as a function of momentum should not be interpreted 
as a dispersion relation of a kind of exciton, because the peaks with broad widths 
constitute an energy continuum generated by a band-to-band transition. We have 
shown that the weight of the — states of the local minority spin in the occupied 
band states strongly depends on the momentumillj. The momentum dependence of 
the weight of the — states of the local minority spin in the band states and 
the dispersion relation of the band states mainly determine the momentum transfer 
dependence of the RIXS intensity in the cuprate and nickelate. 

Although experimental data have been accumulated for the doped cuprates and 
nickelates [2n[22l[23l [^. theoretical analyses are limited on a one-band Hubbard model 
within the exact diagonalization methodplj, and on a three-band Hubbard model 
analysis within the HFA on the basis of the present formalism [25]. An analysis with a 
detailed model like the present paper may be necessary to clarify the momentum and 
incident-photon dependence of the spectra. Since electron correlations are expected 
to be more important in the doped cuprates and nickelates, such studies seem rather 
hard and are left in future. 
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